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ABSTRACT 

We analyse the dependence of the luminosity function of galaxies in groups (LF) on group 
dynamical state. We use the Gaussianity of the velocity distribution of galaxy members as a 
measurement of the dynamical equilibrium of groups identified in the SDSS Data Release 7 
by Zandivarez & Martinez 201 1. We apply the Anderson-Darling goodness-of-fit test to dis- 
tinguish between groups according to whether they have Gaussian or Non-Gaussian velocity 
distributions, i.e., whether they are relaxed or not. For these two subsamples, we compute 
the °- 1 r— band LF as a function of group virial mass and group total luminosity. For massive 
groups, M > 5 x 10 13 M h~ l , we find statistically significant differences between the 
LF of the two subsamples: the LF of groups that have Gaussian velocity distributions have a 
brighter characteristic absolute magnitude (~ 0.3 mag) and a steeper faint end slope (~ 0.25). 
We detect a similar effect when comparing the LF of bright {Mo r ° up — 51og(/i) < —23.5) 
Gaussian and Non-Gaussian groups. Our results indicate that, for massive/luminous groups, 
the dynamical state of the system is directly related with the luminosity of its galaxy members. 
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1 INTRODUCTION 

Group environment plays a key role in the evolution of the 
overall galaxy population in the Universe. Many works in the 
literature have studied galaxies in groups to understand how 
this particular environment affects galaxies and their properties 
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(e.g. iMartmez et all 120021: lEke et all l2004bl: iBalogh et al 



2006; Weinmann et al.l 120061 : lYang et al 



2004 



2009 



Robotham et alj2010h . The understanding of how the different pro 
cesses act upon galaxies requires not only a characterisation of how 
galaxy properties are related to the environment but also to their 
motion therein. The action of the different physical mec hanisms de- 
pend on the dynamics of galaxies within systems (e.g.|Yepes et al.l 
1 199 it iFusco-Femiano & Mencill 19981 : lAbadi et al.ll 19991) . Therels 
evidence that some prop erties of galaxy systems are closely related 
to galaxy dynamics (e.g. | Withmore et al.l 1993; Adam ijit al Jl 19981 ; 



iBiviano et alJl2002l ; lLares et alJl2004l ; lRibeiro et alj20ld) . 

A possible way to characterise the dynamical state of a galaxy 
group is analysing its velocity distribution. It is known that a 
Gaussian velocity distribution is indicative of a group in dynam- 
ical equilibrium, while departures from Gaussianity may indicate 
that perturbative proce sses are working dMenci & Fusco-Femiand 
1996: HoueTai] 2009). However, there is a difficulty in determin- 
ing whether a given velocity distribution differs significantly from 
Gaussian, mainly when studying smaller systems as galaxy groups 
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with o nly a few galaxy members. In a recent work, iHou et all 
(2009) have demonstrated that a reliable distinction can be made 
between Gaussian and non-Gaussian groups even for those with 
low group membership. They conclude that the Anderson-Darling 
(A-D) goodness-of-fit test is the most reliable statistics to distin- 
guish between relaxed and dynamically disturbed systems even for 
those with at least 5 galaxy members. Therefore, this statistical 
method is a very suitable tool to analyse the internal dynamics of a 
system. 

Recently, IZandivarez & Martinez! J201lh (hereafter ZM11) 
used the Seventh Dat a Release of the Sloa n Digital Sky Survey 
(hereafter SDSS DR7; lAbazaiian et al.|[2009h to identify groups of 
galaxies and study several dependencies of the LF. They found that 
the characteristic magnitude brightens and the faint end slope be- 
comes steeper as a function of mass. This change in the luminosity 
function is mainly due to the red spheroids and the varying number 
contributions of the different galaxy types. They also found evi- 
dence of luminosity segregation for massive groups. Moreover, the 
mass trend of the LF is much more pronounced for groups located 
in low density regions. However, the effects of the internal dynam- 
ics of groups on the galaxy LF is an issue that has not been fully 
addressed. Therefore, in this paper we extend the work by ZM1 1 by 
studying the link between the LF and the dynamical state of groups 
by means of their galaxy member velocity distributions using the 
A-D test. The layout of this paper is as follows. In section 2 we 
describe the group sample. The analysis of the LFs is in section 3. 
Finally, in section 4 we discuss the results. 
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Figure 1. The normalised line-of-sight velocity distribution of galaxies in 
Gaussian (upper panel) and Non-Gaussian (lower panel) groups. Solid lines 
are the fitting functions describing the distributions (see text for details). 



2 THE SAMPLE 

For this work, we use the sample of groups constructed by 
ZM1 1. Thi s sample has been i dentified in the Main Galaxy Sam- 
ple (MGS; IStrauss et al.l l2002h of SDSS DR7 which comprises 
galaxies down to an apparent magnitude limit of 17.77 in the r 
ban d. In ZM11, the group ident ification was performed follow- 
ing iMerchan & Zandivare j d2005h : firstly, a standard Friends-of- 
Friends (fof) algorithm links MGS galaxies into groups; and sec- 
ondly, an improvement of the rich group identification is performed 
by means of a second identification on galaxy groups which have at 
least ten members using a higher density contrast. The latter is done 
in order to spl it merged system s or to eliminate spurious member 
detection (see lDfaz et al]|20"05h . The method for estimating group 
centre positions was refined for groups wi th at least ten me mbers 
using an iterative procedure developed bv lDiazetal] d2005h . Due 
to the well known incompleteness of the MGS for r < 14.5, ZM1 1 
excluded galaxies brighter than 14.5. The linking parameters for 
the fof algorithm were set to have a transverse linking length which 
corresponds to a contour over-density of Sp /p = 200 and a line-of- 
sight l inking length of 200 km s _1 . As in lMerchan & Zandivarez! 
(2002), group virial masses were computed as M = a 2 R vlT /G, 
where i? v ir is the virial radius of the system, and a is the veloc- 
ity dispersion of member galaxies ( Limber & Mathews 1960). The 
velocity dispersion a was estimated using the line-of-sight veloc- 
ity dispersion a v , a = y3o\,. To compute a v we used the meth- 
ods described bv lBeersetai] ( ll990h . applying the biweight estima- 
tor for groups with richness N ga i > 15 and the gapper estima- 
tor for poorer systems. The final group sample comprises 15,961 
groups with at least 4 members, adding up to 103,342 galaxies. 
The group sample has a median redshift, velocity dispersion, virial 
mass and virial radius of 0.09, 193 km s"\ 2.1 x 10 13 M Q h~ x , 
and 0.9 h~ 1 Mpc, respectively. 

Galaxy magnitudes used throughout this paper are Pet- 
rosian, are in the AB system and hav e been corre c ted fo r Galac- 
tic extinction using the maps by ISchlegel et alj Jl998t> . Abso- 
lute magnitudes have been computed assuming a flat cosmo- 
logical model with parameters Qq — 0.3, Q,a = 0.7 and 
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Figure 2. The group total absolute magnitude vs. virial mass for ZM11 
groups which have at least 5 members. The thick solid line is the least square 
linear fit between M S ™" P and \og(JVl vir ). Vertical line, M vir = 5 X 
10 13 Mq/i -1 , is the high mass cut-off, while horizontal line, Mf T °' ap — 
51og(h) = —23.5, is the corresponding high luminosity cut-off obtained 
from the estimated linear relation and the high mass cut-off value. 



H = 1 00 h kms^ Mpc" 1 and JsT-corrected using the 
method of iBlanton et alJ {2003a) (KCORRECT version 4.1). We 
hav e also included evolut ion corrections to this magnitude follow- 
ing |Blanton|eral] d2003bh . We have adopted a band shift to a red- 
shift 0.1 for the r band (hereafter 01 r), i.e. to approximately the 
mean redshift of the main galaxy sample of SDSS. 

To distinguish between relaxed groups with Gaussian velocity 
distributions and groups with non-Gaussian dyn amics we adopted 
the A-D goodness-of-fit test. In a recent study, iHou et ai] d2009l) 
have demonstrated that the A-D is one of the most reliable and pow- 
erful tests to measure departures from an underlying Gaussian dis- 
tribution. This test does not require binning or graphical analysis of 
the data. From the outcome of the A-D test, we classify groups into 
two subsamples: the non-Gaussian (NG) groups are those which 
have a confidence level above 90% of not having a Gaussian veloc- 
ity distribution, while Gaussian (G) groups are those whose confi- 
dence level of not having a Gaussian velocity distribution is below 
50%. S ince the A-D te st is reliable when data sets have at least 5 
points dllou et "aT] |2009h . for the purposes of this work we restrict 
the ZM1 1 sample to groups with at least 5 members. From the total 
of 9,387 groups with at least 5 galaxy members, 479 are classified 
as NG while 5,250 as G groups. In Fig.Q]we show, for our samples 
of G and NG groups, the stacked distribution of the radial veloc- 
ity of the galaxies (V) relative to their parent group radial velocity 
(Vcm) normalised to the group velocity dispersion (cr v ). We per- 
form a fitting procedure to clearly describe the velocity distribution 
behaviours. It can be seen (solid lines, Fig. [TJ that the stacked ve- 
locity distribution of G groups is well represented by a Gaussian 
function (upper panel), while the NG groups show clear departures 
from a single Gaussian function, being well fit by the sum of two 
Gaussian functions (lower panel). The small asymmetry in the ve- 
locity profile of NG groups is due to the only 4 groups which have 
more than 90 members. Among them, 3 have left-skewed radial ve- 
locity distribution, and 1 group has a right-skewed one. By exclud- 
ing these large groups, the stacked radial velocity distribution of 
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Figure 3. STY best fitting Schechter function parameters in the r band 
as a function of group mass {left panels) and group total absolute magnitude 
(right panels) for the subsamples of G {filled circles) and NG {open circles) 
groups. Upper panels show the characteristic absolute magnitude, while 
lower panel show the variation of the faint end slope. Vertical error bars are 
the projection of lcr joint error ellipse onto the a and M* axes. Horizontal 
error bars are the 25th and 75th percentiles in each mass/luminosity bin. 

NG groups becomes very close to symmetrical, still non-Gaussian 
and well fit by the sum of two Gaussian functions displaced from 
each other. 



3 THE LF OF GALAXIES IN GROUPS 

This study is based on the analysis of the luminosity function of 
galaxies in groups as a function of a given group physical prop- 
erty, using the group subsamples defined in the previous section. 
As in ZM11, the system physical property adopted is the group 
virial mass. For groups that are not in dynamical equilibrium, i.e. 
those classified as NG, the virial mass might not be a suitable 
measure of the system mass. Thus, comparing the LF of G and 
NG groups as a function of mass can be thought as inappropriate. 
Thus, to complement the analysis of the LF we use also another 
group property which is known to be correlated with mass, the 
group total luminosit y (e.g. iGirardi et al.l2000l : |Popesso et alj200 £ 
iDfaz & Muriell2005l). We com pute group total absolute magnitudes 
following lMoore et alj dl993 ") but using the mass dependence of the 
LF of galaxies in groups as computed by ZM11. We show in Fig. 
[2] the group total absolute magnitudes of the ZM11 groups as a 
function of their virial masses. It can be seen that there is a real 
correlation among these parameters. 

Similarly to the obtained by ZM11 , in all cases, we find that 
the Schechter parametrisation of the LF dSchechtedl 19761) is appro- 
priate for describing the binned LlQ- Therefore, our findings be- 

The binned LF w ere computed using the C~ method jLvnden-Belll97ll ; 
ICholoniews3l987t) . 
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Figure 4. Groups in the single high mass/luminosity bins. Left panels cor- 
respond to groups with M vir > 5 X 10 13 Mq/i — 1 and right panels to 
groups brighter than M® r ° uv — 51og(fc) = —23.5. From top to bottom 
we show the normalised distributions of: virial mass/total luminosity, num- 
ber of members and redshift, of G (solid line) and NG (dotted line) groups. 
Error bars in each histogram are Poisson errors. 

low are expressed in terms of the values of the Schechter function 
shape parame ters, a and M* , on ly, which we compute using the 
STY method dSandage et aO 19791) . 

In Fig. [3] we show the best fitting a and M* parameters of the 
V— band LFs of galaxies in G and NG groups as a function of 
their virial masses (left panels) and total absolute magnitude (right 
panels). We use as many bins as possible to probe the mass/absolute 
magnitude range while having enough data points in each bin to 
produce reliable estimations of the LF. We use 5 and 3 bins which 
include the same number of groups for the G and the NG group 

samples respectively. 

I n agreement with pr evious results (e.g. IZandivarez et al.l 

2006; Robo tham et al.l l2010; ZM11) there is a clear brightening in 
the characteristic magnitude and a decreasing faint end slope as a 
function of mass for both groups subsamples (left panels). As ex- 
pected, the same behaviours are seen as a function of group total 
absolute magnitude (right panels). Comparing the LF parameters 
for G and NG groups, we observe that, within errors, there are no 
significant differences in the Schechter parameters, with the excep- 
tion of the high mass/luminosity tails, where there is an indication 
that galaxies in G groups may have brighter M* and steeper a val- 
ues. 

In order to explore in detail these possible differences at the 
high mass/luminosity tails and to assess their reliability, we re- 
compute the LF for massive/luminous groups by using single bins 
in mass and luminosity. In the single high mass bin we decide 
to include all groups more massive than 5 x 10 13 Mq/i _1 . With 
this choice, we include almost completely the two highest mass 
bins of G groups in Fig. [3] while keeping a number of NG groups 
large enough to have a reliable statistics. For the single high lu- 
minosity bin, we estimate the corresponding value using the previ- 
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Figure 5. STY best fitting Schechter parameters of the 0,1 r band LFs of G (solid line) and NG (dotted line) groups in the single high mass (left panel) and 
the single high luminosity (right panel) bins. We also show their 1, 2 and 3cr confidence ellipses. 



ous high mass cut-off and the least square linear relation between 
group virial mass (log(M v i r )) and group total absolute magnitude 
(Mo, lr — 51og(/i)) for the ZM11 group sample shown in Fig. [2] 
We obtain a group total absolute magnitude cut-off of —23.5 in the 
0,1 r— band. As in the virial mass case, this high luminosity subsam- 
ple includes the two highest luminosity bins of G groups in Fig. [3] 
Now, since we are interesting to demonstrate that the dynamical 
state of the systems is an indicator of the different luminosity be- 
haviours in the high tails, we matched G and NG group distribu- 
tions of virial mass and group total absolute magnitude in each 
subsample (Kolmogorov-Smirnov coefficient > 99.9% in both 
distributions). After performing this procedure we obtain 1225 G 
and 123 NG groups for the high mass subsample and 1373 G and 
197 NG groups for the high luminosity subsample. In the upper 
panels of Fig. [4] we show the distributions of mass (left panel) and 
group total absolute magnitude (right panel) for G (solid line) and 
NG (dashed line) groups in these high mass/luminosity subsam- 
ples. We also show in this Figure the number of members (mid- 
dle panels) and redshift distributions (bottom panels) of G and NG 
groups for both subsamples. It can be seen that the NG groups tend 
have a tail of high membership groups in comparison with the G 
group subsample. This tail can be directly associated with the low 
redshift tail observed for NG groups. 

In the left panel of Fig. [5] we show the STY best fitting 
Schechter parameters for the high mass subsamples of groups along 
with their 1, 2 and 3 a confidence levels. Clearly, the LF of galaxies 
in massive G and NG groups differ at 3a level. The right panel of 
Fig. [5] shows the LF parameters corresponding to G and NG high 
luminosity groups. Again, we observe a clear difference in the LF 
of galaxies in G and NG groups, in this case at a 2a significance 
level. 



4 DISCUSSION 

ZM11 showed that the LF depends not only on local environ- 
ment (group mass, group-centric distance) but also on the large 
scale environment surrounding the groups. Using the same sam- 
ple of groups, in this work we present evidence that the dynami- 
cal state of the system is another important ingredient in the evo- 
lution of the luminosity of galaxies. Our results indicate that, for 
high mass/luminosity groups, the LF of galaxies in G groups have 
brighter M* and steeper a, than the LF of galaxies in NG groups. 
Therefore, the different internal dynamical state of a system is a 
clear indicator of a different history in the galaxy luminosity evo- 
lution. 

Systems of galaxies have Gaussian velocity distributions only 
if they are in dynamical equilibrium. Galaxies in these systems have 
had enough time to suffer the long term action of several physi- 
cal processes during their evolution. Galaxy mergers play a cen- 
tral role in galax y evolution in system s. Some other pr ocesses such 
as str angulation dLarson et alj|l980[). ram pressure |Gunn & Gottl 
1 1972b and galaxy harassment dFarouki & Shapirdll98ll) ~ are more 
efficient in high mass systems, where the effect observed for dif- 
ferent dynamical state of the systems it has been shown to be more 
important. The action of all these processes over the group lifetime 
can produce both, bright and faint galaxies thus providing a plausi- 
ble explanation for our results. 

On the other hand, since the stacked velocity distributions 
for NG groups is well described by two Gaussian functions, it is 
likely that the non Gaussianity is caused by the presence of a mul- 
timodal galaxy population. This behaviour opens the possibility 
for the non Gaussian velocity distributio ns to be the consequence 
of an undergoing merging process (e.g. iMenci & Fusco-Femianol 
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1996) or even multiple merging events (e.g. iGirardi et alj|2005h . 
The different merging populations inhabiting these systems could 
still be experiencing the influence of their own parent halo, and 
hence preserving the galaxy properties corresponding to the indi- 
vidual (smaller) halos that are infalling to form the (larger) non- 
Gaussian group. Therefore, these smaller entities should supply the 
non-Gaussian system with less bright galaxy luminosities that cor- 
respond to less massive/luminous systems. This scenario supports 
the observed fainter characteristic absolute magnitude and the shal- 
lower faint end slope for non-Gaussian systems. Galaxies inhab- 
iting these non relaxed systems are unlikely to feel the influence 
of the environmental physical mechanisms described in the previ- 
ous paragraph, thus preventing the formation of very bright galax- 
ies as we ll as a large number of faint ones. In agreement with this 
scenario, iRibeiro et all J2010h using t he A-D test ov er a sample of 
groups from the 2PIGG catalogue feke et alll2004ah demonstrated 
that galaxies in Gaussian groups are significantly more evolved 
than galaxies in no n Gaussian systems. Also, using a subsample of 
the 2PIGG groups. iRibeiro et all d20 1 lh have shown that non Gaus- 
sian systems are composed of multiple velocity modes, in concor- 
dance with the scenario of secondary infall of clumps at a stage be- 
fore virialisation. Both previous studies were performed analysing 
the surroundings of groups out to 4 times the corresponding radius 
for an overdensity of 200 (47?2oo)- In our work we show that a sim- 
ilar behaviour can be observed from the analysis of the internal dy- 
namics of groups (only galaxy members, mostly inside the virial ra- 
dius) and these different dynamical environments can be evidenced 
in the galaxy luminosities of high mass/luminous systems. Our re- 
sults suggest another way to test models of galaxy evolution, since 
the connection between galaxy luminosities (i.e. astrophysics) and 
the dynamics of the systems should be present. 
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